Optical vortices induced in nonlinear multi-level atomic vapors 
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In a numerical investigation, we demonstrate the existence and curious evolution of vortices in a ladder-type 
three-level nonlinear atomic vapor with linear, cubic, and quintic susceptibilities considered simultaneously 
with the dressing effect. We find that the number of beads and topological charge of the incident beam, as well 
as its size, greatly affect the formation and evolution of vortices. To determine the number of induced vortices 
and the corresponding rotation direction, we give common rules associated with the initial conditions coming 
from various incident beams. © 2012 Optical Society of America 
OCIS codes: 190.4420, 350.2660, 190.6135. 



An optical vortex is an interesting structure that pos- 
sesses a phase defect at a point (called the vortex core) 
and a rotational energy flow around it. In recent decades, 
optical vortices and vortex solitons attracted a lot of at- 
tention from research groups all over the world, for their 
potential applications in optical data storage [1], dis- 
tribution [2] , and processing [3] as well as in the study of 
optical tweezers [4], trapping and guiding of cold atoms 
[5], and entanglement states of photons [6]. To date, re- 
search on optical vortices in many kinds of media, such 
as bulk nonlinear media [7], discrete systems [8], atomic 
vapors [9], dissipative optical systems [10], and Bose- 
Einstein condensates [11], has been reported. In Ref. [12], 
the authors discussed the properties of multidimensional 
beams via atomic coherence. Still, there are many inter- 
esting topics to be researched. 

In this Letter, we consider paraxial propagation 
of probe vortex beams in a ladder-type three-level 
atomic system formed by 35i/ 2 , 3Pi/2 an d 5 Si/ 2 lev- 
els of sodium. The model is described by the nonlinear 
Schrodinger equation (NLSE) of the form [10,12]: 
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where = d xx + d yy is the transverse Laplacian, k 
is the wave number, ip is the amplitude of the beam, 
f3 is the diffusion coefficient, and x is the total sus- 
ceptibility of the atomic vapor system. We consider 
a cubic-quintic NLSE as the underlying propagation 
model that adequately describes the physics of atomic 
vapor systems. The total susceptibility can be obtained 
through the Liouville pathways using perturbation the- 
ory [13]: X = v/K - (v/K^QG, | 2 M + \G 2 \ 2 /d 2 ) 
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is the Rabi frequency of the probe field, G 2 is the Rabi 
frequency of the coupled field, N is the atomic density, 
and Ai ? 2 is the detuning of the probe (coupled) field. 
Tij denotes the population decay rate between the cor- 
responding energy levels \i) and and /iio is the elec- 
tric dipole moment. The second and the third terms in 
X represent the cubic and the quintic contributions to 
the total susceptibility. Equation (1) is similar to the 
complex Ginzburg-Landau equation, with the diffusion 
coefficient j3 coming from the models of laser cavities. 
We launch an incident beam in Eq. (1) of the form (in 
polar coordinates): 



<ip (z = 0, r, 0) =Asech [(r - R ) /r ] 
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with 5 = 0, n 7^ for a necklace, 7^ 0, n 7^ for 
an azimuthon [14], and B = 1 or n = for a vortex. 
Here Ro is the mean radius, ro is the width, A is the 
amplitude, B is the modulation coefficient (1 — B is the 
modulation depth), I is the input topological charge, and 
2n is the number of necklace beads. The beam in Eq. (2) 
can be viewed as a superposition of two vortices with net 
topological charges (NTCs) [15] I + n and I — n, respec- 
tively. The vortex with a larger magnitude of NTC will 
dissipate faster during propagation, due to the diffusion 
term in Eq. (1) [16]. Thus, the vortex with a smaller 
NTC will remain stable for longer during propagation, 
and the overall NTC will correspond to that charge. 

First of all, we set the parameters to be TV = 
10 13 cm" 3 , /iio = 3 x 10~ 29 Cm, A x = 1 MHz, 
A 2 = -1 MHz, G 2 = 40 MHz, r i0 = 2tt x 4.86 MHz, 
r 20 = 2tt x 0.485 MHz, A 600 nm, r 100 //m, 
Ro = 200 /im, and f3 = 0.5 in our numerics for con- 
venience without special statement [17]. The evolution 
of a vortex with I = is shown in the insets in Fig. 
1(a). The input notch of the vortex soon disappears and 
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the vortex changes into a super- Gaussian-like pulse dur- 
ing propagation, with the width and amplitude growing 
dramatically. For comparison, if we set I = 1 and redo 
the evolution, the width and the amplitude still grow, 
but the notch at the origin remains, as shown in the in- 
sets in Fig. 1(b). Even though the amplitude increases 
during propagation, there exists a saturable maximum 
at ~ 10 V/m, as seen from the radial intensity profiles 
in Figs. 1(a) and (b); these results are quite similar to 
the results from Refs. [7] and [12]. The saturation phe- 
nomenon can be explained by the competition between 
the cubic and the quintic nonlinearity [18]. The nonlin- 
earity is defocusing, and as such can support stable vor- 
tices [7] . The beam spreads during propagation, because 
the nonlinearity is too weak to balance diffraction and 
form a soliton. If I is set to 3 and 6, there will be 3 and 6 
notches around the origin, as shown in Figs. 1(c) and (d). 
The right panels there present the corresponding phases 
that demonstrate every notch is a vortex. 




Fig. 1. (color online) (a) and (b) Evolution of vortex 
incidences, for several propagation distances and for I = 

and 1 = 1, respectively, (c) and (d) Output intensities 
(left panels) and phases (right panels) of vortices with 

1 = 3 and I = 6, respectively. The colorbars and physical 
scales of the transverse plane shown here are the same 
in all other figures. 

If we set n = 0.5, we can investigate the evolution of 
a unique "necklace" with crescent shape. In Fig. 2 we 
exhibit a series of evolution snapshots corresponding to 
different I. By comparing Figs. 2(b) with (e), and (c) 
with (d), one can conclude that the number of notches 
that appear in the beam is determined by min{|/ + n|}, 
which is also the absolute value of the NTC of the sur- 
vived vortex component. Even though Figs. 2(a) and (f) 
have the same number of vortices, they rotate in oppo- 
site directions, because the corresponding NTCs are —5 
and 5, respectively. The common rule on how to calcu- 
late the number of vortices and how to determine their 
rotation directions is presented in Table 1, where O and 
O represent rotation senses of the vortices. 

Numerical experiments indicate that there are two 
critical values r cr \ ^42.5 fim and r cr 2 ~ 36.1 fim for the 
beam width ro . The evolution of quadrupole azimuthons 
with n = 2 and different I corresponding to ro > r cr \ 
are shown in Fig. 3. Following the rule from above, the 
number of notches at the origin is determined by NTC; 
so, from (a) to (f) in Fig. 3, the number of notches is 4, 
0, 1, 1, 0, and 4, respectively. However, in (d)-(f), there 
are additional 4 notches around the origin, which equal 




Fig. 2. (color online) (a)-(f) Evolution outputs at z = 
15 cm, from a crescent input (n = 0.5) corresponding to 
I = —5.5, — 1.5, — 0.5, 0.5, 1.5, and 5.5, respectively. 
The top and bottom insets present the input and output 
phases, respectively. 



Table 1. Properties of the induced vortices from necklace 
and azimuthon incidences with different I and n. 



(1) Necklace inputs 


I < — n I = -n -n < I < I 


= 


<l <n 


I = n 


I > n 


O - O 




O 




O 


No. -n-l n + 1 





n — I 





-n + l 


(2) Azimuthon inputs when ro > r cr \ 
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(3) Azimuthon inputs when r < r cr2 with I > n 
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the number of beads. The explanation is that there is 
energy flow around the phase singularities, and at the 
same time the beads fuse due to the diffusion term /3; 
therefore when the speed of energy flow is greater than 
the fusion speed of beads, at the twist of the beads new 
vortices form. It is worth mentioning that the vortices 
around the origin appear as vortex pairs (with charges 
of +1 and —1), so the NTC of the beam is still conserved. 
The common rule for this case is also shown in Table 1, 
under the "Azimuthon inputs when ro > r cr \ case. 




Fig. 3. (color online) Evolutions of azimuthon incidences 
with n = 2, and B = 0.5 for I = -6 (a), -2 (b), -1 (c), 
1 (d), 2 (e), and 6 (f), respectively. The figure setup is 
as in Fig. 2. 

A hexapole azimuthon with I > n is presented in Figs. 
4(a)- (c); its ro = 40 fim is chosen to fall in-between r cr 2 
and r cr \. According to the rule from Table 1, the to- 
tal number of induced vortices should equal n + l = 12. 
However, from Fig. 4 we see that there are more and 
more vortices appearing during propagation. The reason 
is that the speed of energy flow is the largest, as com- 
pared to those of fusion and spreading. From the phases 
at different propagation distances, we see that the en- 
ergy flow brought by the vortices around the origin of 
the fused beam is always faster than that at the edge 
of the beam; therefore, the asynchrony preferably forms 
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new phase singularities at the edge of the beam, and 
new vortices are induced correspondingly. But we can- 
not give a certain rule for this case, because the number 
of induced vortices is greatly affected by the beam width. 
What one sees in Fig. 4(c) is 6 vortices at the core plus 
12 pairs of ±1 vortices induced at the rim of the beam. 

If we set ro = 30 /im, which fulfills the condition 
ro < r cr 2, and redo the propagation of the beam used 
in Figs. 4(a)- (c), we find that the number of induced 
vortices is 6, which can be calculated from I — n = 6, 
as shown by the output intensities and phases in Fig. 
4(d). The common rule for this case is exhibited in Ta- 
ble 1 under "Azimuthon inputs when ro < r cr 2 with 
I > n" (other cases are the same as those under "Az- 
imuthon inputs when ro > r cr i"). The lower the number 
of beads, the more energy in each of the beads, which 
will strengthen the energy flow. That is why the number 
of induced vortices is I + n if n = 1. However, the exact 
rule for n = 2 is not certain, because the energy flow is 
weakened, and whether it can produce more vortices or 
not, it is up to the value of I set for the initial beam. If I is 
bigger, the fusion will be accelerated, so the production 
of new vortices will be limited. Here, I = 5 is a boundary 
for this case: the rule is I — n for I > 5 and I + n for I < 5. 
Figures 4(e) and (f) display two numerical experiments 
corresponding to I = 9 and I = 4, in which the number 
of vortices are I — n = 7 and I + n = 6, respectively. 




Fig. 4. (color online) (a)-(c) Evolution of an azimuthon 
with I = 9, n = 3, B = 0.5, ro = 40 /im at differ- 
ent propagation distances, (d)-(f) show azimuthons with 
ro = 30 /im at 30 cm with (I = 9, n = 3), {I = 9, n = 2), 
and (I = 4, n = 2), respectively. The top row are the 
intensities, the bottom row are the phases. 

In conclusion, we have demonstrated that optical vor- 
tices can form from vortex, necklace, and azimuthon inci- 
dences with different topological charges, in multi-level 
atomic vapors when linear, cubic and quint ic suscepti- 
bilities are considered simultaneously. The appearance 
of vortices results from a combined action of the num- 
ber of topological charges, the beam width of incidences, 
the diffusion effect, cubic-quintic nonlinearities, and the 
loss/gain in the medium, simultaneously. We have for- 
mulated common rules of finding the number as well as 
the rotation direction of the induced vortices. The effects 
can be observed in sodium as well as other atomic va- 
pors with incident beams produced by using multi-wave 
interference [19] or phase mask [20] methods in a similar 
experimental setup as that in Ref. [21]. Our findings open 



a new venue and introduce a new method for studying 
vortices, and may broaden the field of their applications. 
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